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Low-dimensional molybdenum and tungsten oxides 
and bronzes have been the focus of much attention be- 
cause of the charge density wave (CDW) instabilities they 
exhibit [2| . An important step in our present under- 
standing of the physics of low-dimensional metals was 
the discovery of non linear transport due to sliding of 
the CDW in the K . 3 MoO 3 blue bronze 3]. This re- 
port launched a wide interest on these materials and, 
as a result, a number of quite intriguing discoveries re- 
lated to the condensation of CDWs were reported p], Q . 
The CDW periodicity is incommensurate with the lattice 
constant in the blue bronzes or the 77— and 7— M04O11 
Magneli phases but commensurate in the sodium and 
potassium purple bronzes. Another widely studied family 
of these compounds is that of the monophosphate tung- 
sten bronzes, which exhibit successive CDWs with mod- 
ulations varying considerably for the different members 
of the family UQ. The condensation of these CDWs is 
usually driven by instabilities of the Fermi surface (FS) 
[f| . In some cases the FS clearly exhibits nesting features 
which are at the origin of the CDW instability. However, 
in other cases the FSs are apparently non nested and for 
some time it was difficult to understand the origin of the 
CDWs. The development of the so called hidden nesting 
concept [H provided a simple way for analyzing the ap- 
parently non nested FSs in terms of a superposition of 
nested ones. 

Recently, Raj et al. 0] have reported angle- 
resolved photoemission spectroscopy (ARPES) and first- 
principles density functional theory (DFT) results con- 
cerning the FS of the hexagonal potassium tungsten 
bronze, K0.25WO3 [8] (see Fig. QJa)). On the basis of 
their results they have suggested that the well known 
high-temperature resistivity anomaly of this phase origi- 
nates from a possible CDW instability as a result of the 
existence of hidden pseudo one-dimensional (ID) bands. 
However, this conclusion is incorrect since it is based on 
the analysis of data for the T-K-M section only of the 
hexagonal Brillouin zone (see Fig. H^b)). Here we re- 
port first-principles DFT results showing that K0.25WO3 
does not possess pseudo- ID partially filled bands and is 
instead a 3D metal. In addition, we analyze how the 
crystal structure determines the nature and shape of the 
partially filled bands and thus the transport properties. 

The present calculations were carried out using a nu- 
merical atomic orbitals DFT Q approach implemented 
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FIG. 1: (Colour online) (a) Crystal structure of the K0.25WO3 
system; (b) Hexagonal Brillouin zone of K0.25WO3; (c) Cal- 
culated band structure for K0.25WO3 where the dashed line 
refers to the Fermi level. 



m the SIESTA code pjj. We have used the generalized 
gradient approximation to the exchange and correlation 
potential and, in particular, the functional of Wu and 
Cohen Only the valence electrons are considered 

in the calculation, with the core being replaced by norm- 
conserving scalar relativistic pseudopotentials [12| factor- 
ized in the Kleinman-Bylander form pjjj]. We have used 
a split-valence double-^ basis set including polarization 
orbitals for all atoms, as obtained with an energy shift of 
0.02 Ry [H. 

The energy cutoff of the real space integra- 
tion mesh was 300 Ry and the Brillouin zone (BZ) was 
sampled using a grid of (12x12x12) fc-points [15|. The 
existence of partially empty potassium sites was taken 
into account via the virtual crystal approximation. 

The essence of the reasoning by Raj et al. [7] is that 
both the calculated and ARPES (a* 6*) sections of the 
FS result from the weak hybridization of pairs of parallel 
lines, perpendicular to the (a* +6*) direction and those 
equivalent by 120° rotations. From this observation they 
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conclude that "the FS is the consequence of hidden (ID) 
one-dimensional bands". This clearly means that these 
bands should have a nil dispersion along the c* direction. 
The calculated band structure for K0.25WO3 in the re- 
gion around the Fermi level is shown in Fig. [Tfc. It is 
immediately evident that the three partially filled bands 
exhibit a very large dispersion along c* (i.e. along T — > 
A, M — > L and K — > H). Consequently, there are no 
partially filled ID bands along the direction (a* +6*) and 
those equivalent by 120° rotations. Since the three par- 
tially filled bands are already slightly above the Fermi 
level at the K point, this means that their dispersion 
in the c*=0.0 section is also large and consequently, the 
partially filled bands of K0.25WO3 should be 2D or 3D. 

The calculated FS for K0.25WO3 is reported in Fig. [2] 
for different sections perpendicular to the c* direction. 
The FS in the plane at c*=0.0 is practically identical 
with the theoretical and experimental sections reported 
by Raj et al. 0] (the very slight differences are due to 
the different way in which the potassium doping has been 
treated in the calculations: via a rigid band approach 
using hexagonal WO3 by Raj et al. 0] or through the 
explicit consideration of K cations via a virtual crystal 
approximation here). This section may indeed be seen 
as resulting from the interaction of three pairs of paral- 
lel lines. However, it is also clear that when departing 
from the plane at c*=0.0 the shape of the FS section se- 
riously changes. In the plane at c*=0.25 any trace of the 
apparent hidden ID bands is already lost and none of 
the sections at c* > 0.3 contributes to the FS. Thus, the 
three components of the FS are 3D as expected from the 
calculated band structure. In fact, this conclusion could 
have been reached by Raj et al. [7] since they noted that 
there is no band crossing the Fermi level in the AHL 
plane. Alternatively, simple electron counting arguments 
also lead to the same conclusion. The K0.25WO3 unit cell 
contains six formula units so that there are 1.5 electrons 
to fill the three t2^-type bands of Fig.[]Jc). Consequently, 
if three bands were ID they should be one-quarter filled 
and the nesting vectors can not be those proposed by Raj 
et al. [7] which would be appropriate for half-filled ID 
bands. 

Coming back to the FS section of Fig. [2j let us note 
that fiat portions with weak curvature may remain until 
sections with c* ~ 0.1. Therefore, it can not be com- 
pletely discarded that weak maxima of the electron-hole 
Lindhard response function can arise. However, even as- 
suming that this weak feature could stabilize a CDW 
modulation [16[, the wave vectors of the modulations 
proposed by Raj et al. [7] are not those which would 
be expected. In fact, when a FS like that of Fig. (2][a) 
results from the weak interaction of ID bands, the ap- 
propriate nesting vectors, i.e. those leading to strong 
maxima in the Lindhard function, are those which allow 
the simultaneous nesting of two different pairs of parallel 
lines. This is the so called hidden nesting [6[ mechanism 
underlying the CDW instabilities exhibited by a num- 
ber of oxides and bronzes like the potassium and sodium 
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FIG. 2: Calculated FS for K0.25WO3 for different sections of 
the Brillouin zone perpendicular to the r — > A direction (i.e. 
c*): (a) 0.0c*, (b) 0.0833c*, (c) 0.1666c* and (d) 0.25c*. 



purple bronzes, the 77— and 7— M04O11 Magneliphases 
and several monophosphate tungsten bronzes p], d, 0, [|| . 
In the present case, the nesting vectors would be l/2a*, 
1/26* and (l/2a*-l/26*). As mentioned above, although 
the condensation of CDW modulations of this type can 
not be completely discarded, the driving force should be 
quite weak. Let us mention that the super lattice spots re- 
ported by Krause et al. [17[ are quite sample dependent. 
Moreover, although the nesting vectors proposed by Raj 
et al. 0] may be found in the diffraction patterns, other 
wave vectors also exist [l7| . In view of these observations 
we believe that there is no compelling reason to believe 
that the superlattice spots and the weak resistivity cusp 
at 350 K [17] are due to the condensation of CDWs. We 
find the arguments of Krause et al. [17[ concerning the 
ordering and migration of potassium vacancies more con- 
vincing. Also, the FS determined by photoemission at 20 
K in the (a* 6*) reciprocal plane [7[ does not reveal the 
formation of partial gaps expected from the stabilization 
of a CDW. 




FIG. 3: (Colour online) (a) Schematic diagram showing that 
every hexagonal octahedral plane can be seen as resulting 
from the intersection of octahedral chains along the a, b and 
(a+6) directions, (b) Top view of a xz/yz orbitals chain along 
the (a-\-b) direction (see text). 
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With this analysis in mind it is puzzling that both 
the experimental and theoretical section of the FS for 
c*=0.0 is made of sets of weakly interacting pairs of par- 
allel lines, as if it was originating from partially filled 
ID bands. This feature is however easy to understand 
after analysis of the band structure The three par- 
tially filled bands are essentially built from the tungsten 
xz and yz orbitals which mix with the oxygen p orbitals 
in an antibonding way. In an ideal hexagonal structure, 
the tungsten xz and yz orbitals are antisymmetric with 
respect to the basal planes of the octahedra. For wave 
functions corresponding to a 0.0 c* component of the k- 
vector, the p XjV orbitals of the apical oxygens can not mix 
by symmetry into the wave function. Consequently, there 
is no interaction between the tungsten xz/yz orbitals of 
different hexagonal planes for this specific section of the 
Brillouin zone. In this particular case, the band disper- 
sion can only arise from interactions within a hexagonal 
plane. The pseudo ID nature of the bands for this sec- 
tion is easy to understand when it is realized that: (a) 
as shown in Fig. [31(a), one hexagonal octahedral plane 
can be seen as resulting from the intersection of octahe- 
dral chains along the a, b and (a+b) directions, in such 
a way that every octahedron belongs simultaneously to 
two different types of chains, and (b) every octahedron 
possesses one pair of xz/yz orbitals which can mix and 
rotate leading to two new orbitals avoiding as much as 
possible the antibonding mixing with the basal oxygen 
p z orbitals, thus leading to the lowest energy t2^-based 
crystal orbitals. In other words, every octahedron can use 
independently the two appropriate xz/yz combinations 
to make extended interactions along the two chains com- 
mon to every octahedron. This leads to chains of orbital 
interactions like that schematically shown in Fig. [31(b) 
(top view), which run along the (a+6) direction. These 
chains of orbitals do not interact with orbitals of parallel 
chains because the tungsten xz/yz combination can not 
overlap with the p z orbital of the oxygen atoms placed 
at the basal site perpendicular to the chain direction as 



a result of the local symmetry (see Fig. [3]d) . Hence, the 
absence of such oxygen contributions isolates different 
parallel chains. In addition, the fact that two types of 
chains intersect at each octahedron does not lead to a 
sizeable interaction because of the local orthogonality of 
the two combinations of xz/yz orbitals. Taking into ac- 
count the actual distortions of the lattice with respect 
to the ideal one does not alter in any significant way 
this description of the partially-filled bands. As a matter 
of fact, the crystal orbitals at T are essentially made of 
sets of parallel chains of orbitals like that of Fig. 02(b) 
(or those equivalent along the a and b directions) which 
practically do not communicate. When the component 
of the fc-vector changes along the direction of the chain, 
the p z orbitals of the oxygens shared along the chain 
can mix into the crystal orbital because of the associated 
phase changes and the band acquires dispersion. How- 
ever, when the k- vector component changes along the in- 
terchain direction, the chains practically do not commu- 
nicate because of the local symmetry argument pointed 
above. This means that for this particular section of 
the Brillouin zone the system can be considered to result 
from the very weak hybridization brought about by the 
octahedral distortions of three ID chains running along 
the a, b and (a+6) directions. However, as soon as the c* 
component of the /c-vector departs from zero the inter- 
octahedral planes interactions are switched on and any 
trace of the pseudo- ID behavior is lost. In summary, we 
find that K0.25WO3 does not exhibit ID partially filled 
bands and that even if the Lindhard function could ex- 
hibit weak maxima due to nesting of the FS, the more 
intense ones should not be those previously proposed Q • 
We believe that arguments based on potassium vacancy 
ordering are more likely when trying to explain the origin 
of the resistivity anomaly of this bronze. 
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